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Abstract 

A collocation method is presented for numerical solution of a typical 
integral equation Rh := f D R(x,y)h(y)dy — f(x), x € D of the class 
TZ, whose kernels are of positive rational functions of arbitrary selfad- 
joint elliptic operators defined in the whole space R r , and D C R r is a 
bounded domain. Several numerical examples are given to demonstrate 
the efficiency and stability of the proposed method. 

MSG: 45A05, 45P05, 46F05, 62M40, 65R20, 74H15 

Key words: integral equations in distributions, signal estimation, collocation 
method. 

1 Introduction 

In [4] a general theory of integral equations of the class 1Z was developed. The 
integral equations of the class 1Z are written in the following form: 

Rh:= f R(x,y)h{y)dy = f(x), xeD:=DUT, (1) 

where D e W is a (bounded) domain with a (smooth) boundary T. Here the 
kernel R(x, y) has the following form (4) [3 El [7] : 

R(x,y) = f P(X)Q- 1 (XMx,y 1 X)dp(\), (2) 

J A 



where P(A), Q(X) > are polynomials, degP = p, degQ — q, q > p, and $, p, A 
are the spectral kernel, spectral measure, and spectrum of a selfadjoint elliptic 
operator C on L 2 (W) of order s. It was also proved in [3] that R : H~ a (D) — ► 
H a {D) is an isomorphism, where H a (D) is the Sobolev space and H~ a {D) its 
dual space with respect to the L 2 (D) inner product, a = s ( q ~ p > , Here the space 
H~ a (D) consists of distributions in H a (W) with support in the closure of D. 
In this paper we consider a particular type of integral equations of the class 
Tl with D = (-1,1), r = 1, C = -id, d := £, A e (-00,00), dp(A) = cZA, 

A) = P(A) = 1, Q(A) - ^±1, 8 = l,p = 0, q = 2 and a = 1, 

i.e., 

Rh(x) := y e-l-*l%)di/ = /(as), (3) 

where /i G 1,1] and / G ff 1 [ — 1,1]. We denote the inner product and 

norm in iP[— 1, 1] by 



(u,v ), := 



fu(x)v(a;) + u'(x)w'(a;)dx^ <fx b,u6 1, 1]), (4) 



and 

H?:= £(KaOI 2 + l«'(*)l 2 )<fc, (5) 

respectively, where the primes denote derivatives and the bar stands for complex 
conjugate. If u and v are real- valued functions in H [— 1, 1] then the bar nota- 
tions given in (0} can be dropped. Note that if / is a complex valued function 
then solving equation Q is equivalent to solving the equations: 

J e-\ x -y\h k (y)dy = f k (x), fc = l,2, (6) 

where h\{x) :=Reh(x), h,2(x) :=lmh(x), fi(x) :=Ref(x), f%(x) :=lmf(x) and 
h(x) — hi(x)+ifi2(x) 1 i = v 7 — 1- Therefore, without loss of generality we assume 
throughout that f(x) is real-valued. 

It was proved in [5] that the operator R defined in ([3]) is an isomorphism 
between £T _1 [— 1, 1] and i? 1 ]— 1, 1]. Therefore, problem ([3]) is well posed in the 
sense that small changes in the data f(x) in the iP[— 1, 1] norm will result in 
small in H~ x [—1, 1] norm changes to the solution h(y). Moreover, the solution 
to (|3]) can be written in the following form: 

h(x) = a-i5(x + 1) + a S(x - 1) + g(x), (7) 

where 
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and 5(x) is the delta function. Here and throughout this paper we assume that 
/ € C a [-l,l], a > 2. It follows from © that h(x) = g(x) if and only if 
/(-l) = /'(-!) and /(1) = -/'(1). 

In [7] the problem of solving equation © numerically have been posed 
and solved. The least squares method was used in these papers. The goal of this 
paper is to develop a version of the collocation method which can be applied 
easily and numerically efficiently. In 18] some basic ideas for using collocation 
method are proposed. In this paper some of these ideas are used and new ideas, 
related to the choice of the basis functions, are introduced. In this paper the 
emphasis is on the development of methodology for solving basic equation (1) 
of the estimation theory by a version of the collocation method. The novelty of 
this version consists in minimization of a discrepancy functional (26), see below. 
This methodology is illustrated by a detailed analysis applied to solving equation 
(3), but it is applicable to general equations of the class 7Z. One of the goals 
of this paper is to demonstrate that collocation method can be successfully 
applied to numerical solution of some integral equations whose solutions are 
distributions, provided that the theoretical analysis gives sufficient information 
about the singular part of the solutions. 

Since / e C Q [-1, 1], a > 2, it follows from © that g € C[-l, 1]. Therefore, 
there exist basis functions fj{x) G C[— 1, 1], j = 1, 2, . . . , m, such that 

max \g(x) — g m {x)\ — > as m — > oo, (10) 

!B£[ — 1,1] 

where 

m 

g m {x) := ^c^pjix), (11) 
i=i 

C 7™' ) J = l,2,...,m, are constants. Hence the approximate solution of equation 
([3]) can be represented by 

h m (x) = £$5(x + 1) + c ( ™ ] 5{x - 1) + g m (x), (12) 

where c*- , j = —1,0, are constants and g m (x) is defined in (fTTj) . The basis 
functions tpj play an important role in our method. It is proved in Section 3 that 
the basis functions tpj in (|12[) can be chosen from the linear B-splines. The usage 
of the linear B-splines reduces the computation time, because computing (|12[) 
at a particular point x requires at most two out of the m basis functions ipj . For 
a more detailed discussion of the family of B-splines we refer to |10) . In Section 
2 we derive a method for obtaining the coefficients Cj , j = —1,0,1, ... ,m, 
given in (|12[) . This method is based on solving a finite-dimensional least squares 
problem ( see equation (f53"f below ) and differs from the usual collocation method 
discussed in [2] and [3]. We approximate — Rh m \\ \ by a quadrature formula. 
The resulting finite-dimensional linear algebraic system depends on the choice 
of the basis functions. Using linear B-splines as the basis functions, we prove 
the existence and uniqueness of the solution to this linear algebraic system for 
all m = m(n) depending on the number n of collocation points used in the left 
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rectangle quadrature rule. The convergence of our collocation method is proved 
in this Section. An example of the choice of the basis functions which yields the 
convergence of our version of the collocation method is given in Section 3. In 
Section 4 we give numerical results of applying our method to several problems 
that discussed in [7]. 



2 The collocation method 



In this Section we derive a collocation method for solving equation @. From 
equation §5§ we get 

Rh(x) = a_ie- (3;+1) + a e- {1 - x) 

r f 1 \ ( 13 ) 

-* J e«g(y)dy + e* j e-y g (y)dyj = f(x). 

Assuming that / € C 2 ([— 1, 1]) and differentiating the above equation, one ob- 
tains 

(Rh)'(x) = -a_ie- (x+1) +a e- (1 - x) 

+ (e x t e-yg{y)dy - e~* f e v y{y)dy\ = f'(x). ' ' ' ' 



Thus, fix) and f'(x) are continuous in the interval [—1,1]. Let us use the 
approximate solution given in (|T2j) . From (fT3| . ([T4|) and (fT2|) we obtain 



Rh m {x) = c^e-^ + c^e-^ 



4 m) [ e- x / e v Vj{y)dy + e x [ e-* Vj (y)dy) := / m (x), ' ' ' ' 

j'=i 



and 



3=1 



„( m ) 



e v ipj(y)dy-e x / e v ipj(y)dy 



(fm)'(x). 

(16) 



Thus, Rh m {x) and (Rh m )' (x) are continuous in the interval [—1, 1]. Since f(x) — 
Rh m {x) and /'(a;) — (Rh m )'(x) are continuous in the interval [—1,1], we may 
assume throughout that the functions 



and 



Ji,m ■= [f(x) - Rh m (x)] 2 
J2, m ■= lf'(x) - (Rh m )'(x)} 2 



(17) 
(18) 
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are Riemann-integrable over the interval [—1,1]. 
Let us define 

q-i(x) :=e-(* +1 >, q Q (x) := e~^~ x \ 
q 3 (x) := J e-\ x -y\tpj{y)dy, .7 = 1,2,..., 

and define a mapping W. n : C 2 [— 1, 1] — > R™ by the formula: 



m, 



(19) 



/ 0(*i) \ 



4>{x) & C 2 [-l,l], 



(20) 



where 



Z2 



■ = < 



e R" | := zfe), z e C 2 [-l,l] 



(21) 



and Xj are some collocation points which will be chosen later. We equip the 
space R™ with the following inner product and norm 



i=i 



(22) 
(23) 



II^m^E^K + K-) 2 ], «eR?, 

respectively, where Uj := u(xj), u'j := u'(xj), Vj := w(a;j), := v'(xj), and 
Wj > are some quadrature weights corresponding to the collocation points Xj , 
j = 1,2, ...,n. 

Applying to i?/i m , one gets 



E- 



-(l/-a:i) 



V>j{y)dy 1 i = 1,2,..., n, 



where m = m(n) is an integer depending on n such that 
m(n) + 2 < n, lim m(n) = oo. 

Let 



(24) 

(25) 
(26) 
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where / and Rh m are defined in ^ and (|T5|) . respectively, H n defined in (|2"U)) , 

/ O-X \ 

CO 
Cl 



c„) i defined in J23J) and c< m ) 



G 



pm+2 



Let us choose 



\ C m J 



(«) 2 • 10 

w- = -, j = 1, 2, . . . , n, 

J n 



(27) 



and 



i,- = -1 + (i - l)s, s:=-, j = l,2,...,n, (28) 
so that ||% n (/ — -R^m)|| 2 (n) ! is the left Riemannian sum of ||/ — Rh m \\\, i.e., 
HI/ - Rh m \\l - \\H n (f - Rh m )\\ 2 w{n) t \ := S n -> as n -> 00. (29) 



Remark 2.1. // Ji !m and J2. m are in C [— 1, 1], where J\. m and J2,m 
defined in (|1T[) and (|18|) . respectively, then one may replace the weights 
with the weights of the compound trapezoidal rule, and get the estimate 

-1 



» 



S n = 



(Jl,m(x) + J 2 ,m{x))dx - ^Wf' (J-]_, m (Xj) + J 2 , m (Xj)) 
3=1 



(30) 



where S n is defined in (I29[) and 

A/ := I J' 1<m (l) + Ji, m (l) ~ W, m (-1) + ( 31 ) 
.Here we /lave used £/ie following estimate of the compound trapezoidal rule J7J 



?7(x)dx - 

3=1 



< 



(b-a) 2 
Yin 2 



r\ (x)dx 



Yin 1 



(32) 

where r\ € C 2 [a,b]. Therefore, if Dj < C for all m, where C > is a constant, 
then 5 n = O (■\) . 



The constants Cj in the approximate solution h m , see ([12|) . are obtained 
by solving the following least squares problem: 



minG n (c< m >), 



(33) 



where G„ is defined in ([26 
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Q = j^ w ^( Em , d ^+ElJ-^A fc = -l,0,l,...,m, (34) 



A necessary condition for the minimum in (1331) is 
>E m j E , dE' m 

where 

E m>l := (f - Rh^ixi), E' m>l :={f -Rh m )'{xi), 1 = 1,2,..., n. (35) 

Necessary condition (|3"4")l yields the following linear algebraic system (LAS): 

A m+2 c {m) = F m+2 , (36) 

where £ R m+2 , A m +2 is a square, symmetric matrix with the following 
entries: 



(A m+2 ) hl :=2j2wl n) e- 2 ^ +1 \ (A m+2 ) h2 = 0, 

n 

(A m+2 )i.j :=2Y J w\ n) Ci d - 2 e-^ l+l \ j = 3,...,m + 2, 
i=i 

n 

{A rn+2 ) 2a :=2^ Wi ( " ) e - 2 ( 1 -^), 

1=1 

n 

(A m+2 ) 2 ,j := 2 wl n) B l j- 2 e-V- xt \ j = 3,...,m + 2, 



-i 



i=i 

i = 3, . . . , m + 2, j = i, . . . , m + 2, 

(A m+2 )j,i = {A m+2 )ij, i,j = 1,2,..., to + 2, 

is a vector in R m+2 with the following elements: 



(JW 2 )i := £>< n) (/(*;) - f'( Xl ))e-^ = (H n Rh,Hq-i) w(n)A 
i=i 



(37) 



(38) 



(39) 



(40) 



1=1 

n 

(F m+2 )i ~J2 w i n) [f(xi)(Ci,i- 2 + Bi,i_ 2 ) + f(xi)(B hi _ 2 - C M _ 2 )] 
1=1 

= (HnRh,Hqi} w(n)1 , i = 3,...,m + 2, 

and 

#U := / e-b-^ipjWdy, C u := / e^'"^^)^. (41) 

J-i 
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Theorem 2.2. Assume that the vectors H n <lj> j = —1,0, 1, . . . , m are linearly 
independent. Then linear algebraic system (|36[) is uniquely solvable for all m, 
where m is an integer depending on n such that (|25p holds. 

Proof. Consider qj G 1, 1] defined in (fH?)) . Using the inner product in R™, 

one gets 

{A m+2 )i,j = (H n q i -2,'r{nqj-2) w (n) tV i, j = 1, 2, . . . , m + 2, (42) 

i.e., A m+ 2 is a Gram matrix. We have assumed that the vectors H n qj G R™, 
j = —1,0,1,.. are linearly independent. Therefore, the determinant of the 
matrix A m+ 2 is nonzero. This implies linear algebraic system (|36[) has a unique 
solution. 

Theorem 12.21 is proved. □ 

It is possible to choose basis functions <pj such that the vectors W n qj, j — 
— 1, 0, 1, . . . , m, are linearly independent. An example of such choice of the basis 
functions is given in Section 3. 

Lemma 2.3. Let y m := c^™) be the unique minimizer for problem (|33[) . Then 



where G n is defined in (|26|) and m is an integer depending on n such that (|25p 



Proof. Let 

h(x) = a_i<5(x + 1) + aoS(x — 1) + g(x) 
be the exact solution to ([3]), Rh = f, where g(x) E C[— 1, 1], and define 



G n (y m ) ->• as n -> oo, 



(43) 



holds. 



hm (^) 



a_i<5(x + 1) + ao<S(a; - 1) + 3 m (x), 



(44) 



where 




(45) 



Choose g m {x) so that 



max lo(x) 



ft: 



'm 



(x)| — > as m — > oo. 



(46) 



Then 

Gn(ym) < \\H n (f - Rh m )\\l 



(47) 



because ?/ m is the unique minimizer of G„ . 

Let us prove that \\% n (f — Rh m )\\ 2 (n) 1 — ► as n — > oo. Let 



: = ~ Rh m {x), 



W 2 , 



m 



f( x ) - (Rh m y(x). 
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Then 



W hm (x) = e- X J ^y(g(y)-~g m (y))dy + e x J e~* (g(y) - ~g m {y))dy (48) 
and 

W 2 , m (x) = e x j e-y(g(y) - ~g m (y))dy - e~* J" e%(y) - g m (y))dy. (49) 

Thus, the functions [Wi jm (a:)] 2 and [W2,m(a0] 2 are Riemann-integrable. There- 
fore, 

6 n := HI/ - - - Rh m )\\ 2 wWil \ ->■ as n -> oo. (50) 

Formula (|50p and the triangle inequality yield 

\\H n (f - Rhm)\\l Mil <S n + \\f- Rh m \\l (51) 

Let us derive an estimate for ||/ — Rh m \\\ . From (14"5|) and (|49l) we obtain the 
estimates: 



\W x , m {x)\ < max ^(y) - g m („)| e~ x / e^dy + / e -»d» 

V J-i Jx 

= rzs f \g{y)-g m {y)\[e-*{e*-e- l ) + e*{e-*-e- l )\ (52) 

= max ^(y) - § m (y)| [(2 - e" 1 - - e" 1 ^)] < S mA 
«e[-i.i] 



and 



|W 2>m (a;)| < y m«K 1] Itffc/) - ffm(l/)l ( / e^dy + e~ x / e^rfy ) < 5 m ,x, 



where 



(53) 



5 m> i:=2 max | g(y) - g m (y) \ . (54) 
»e[-x,i] 



Therefore, it follows from (1521 and (j5"3")l that 



\\f-Rh m \\i = J \W hm {x)\ 2 dx+ I \W 2 , m (x)\'dx 
< 4€ i, 



(55) 



where is defined in ff54|) . Using relation (j46]), we obtain lim m _>.oo <5 m ,i = 0. 
Since m = m(n) and limn-^ m(n) = oo, it follows from (|51[) and ()55[) that 
||^ n (/ ~ ^m)||w,i ^ as 7i ^ oo. This together with flU} imply G n {y m ) -> 
as n — > oo. 

Lemma T2. 31 is proved. □ 
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Theorem 2.4. Let the vector 



system ([55]) and 



( \ 

r ( m ) 

„(">) 



£ K m+2 solve linear algebraic 



h m (x) = c [ ^5{x + 1) + ci m) S(x - 1) +Y / c ( f l) ^(x). 



Then 



\\h — h m \\ H -i — > as n — > oo. 

Proof. We have 

11^ — ^ra II H- 1 [-1,1] = 11-^ (/ — ^m)llij-i[-i,i] 

— 1 II 2 



< \\R 



Iff 1 [— 1,1]— 1,1] II f -R^mUl 

< C (G„(4t) + |||/ - Rh m \\\ - G„(4t)|) 

< C[G n (<^L) + 5 n ] -»■ as n -> oo, 



where C > is a constant and Lemma T2. 31 was used. 
Theorem 12.41 is proved. 



(56) 



(57) 



□ 



3 The choice of collocation points and basis func- 
tions 

In this section we give an example of the collocation points X{, i = 1, 2, . . . , n, 
and basis functions cfj, j = 1,2, ...,m, such that the vectors H n qj, j = 
— 1, 0,1, ... ,m, are linearly independent, where m + 2 < n, 



p— 1— X 2 



1+^2 



-1+K„ 



(58) 



and 



rL n qj 



e -x2 e y ipj(y)dy + e X2 $\ e- y ip j (y)dy 
\ e~ Xn ll 1 e y (p 3 (y)dy + e Xn J^ 1 e-yip 3 (y)dy J 



, j = l,2,...,m. (59) 
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<Pm{x) 



^2j-3 x 2j-2 %2j-l x 2j x 2j+l 



Figure 1: The structure of the basis functions ipj 



Let us choose wj™' and Xj as in (|27)) and (|28|) , respectively, with an even 
number n > 6. As the basis functions in C[— 1, 1] we choose the following linear 
B-splines: 



tpi{x) 
tpj(x) 



ip 2 (x) Xi < X < x 3 , 
0, otherwise, 

- (j - l)2s), %2j-3 < x < X 2 j-1, 
ip2(x - (J - l)2s), X 2 j-1 <X< X 2 j+1, 
0, otherwise, 



(60) 



3 = 2,...,m- 1, 



ipi(x — (m — l)2s), x n -\ < x < 1, 



0, 



otherwise, 



where 



2 +1 ' 



2 

s := -, 
n 



and 



ip 2 (x) := 



x — xi + 2s 
2s" ' 
— (x — xi 



2s) 



(61) 



2s 



Here we have chosen x 2 j-i, j — 1, 2, . . . , ^ + 1, as the knots of the linear IB- 
splines. From Figure 1 we can see that at each j = 2, . . . , m— 1, <Pj(x) is a "hat" 
function. The advantage of using these basis functions is the following: at most 
two basis functions are needed for computing the solution h m {x) 1 because 



km {x) 



„( m ) 



X21-1, 



c\ 'ipi(x) + c\ + {(pi + i{x), X 2l -i < X < X21+1, 
From the structure of the basis functions ipj we have 
<Pi(x) =0, CC3 < x < 1, 

<Pj(x) = 0, —1 < x < x 2 j-3 and x 2 j+\ < x < 1, j = 2,3, 
=0, -1 < x < x„_i. 



Z = 2,. 



4 <«> 



(63) 
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Let ('HnQj)i be the i-th element of the vector Tl n qj, j = — 1, 0, 1, . . . , to. Then 

{H n q-x)i = e-^ +Xi \ i = l,2,...,n, (64) 

CH n qoh = e-^- x '\ i = l,2,...,n. (65) 
Using (j6"3"j) in ([5^]) . we obtain 

f i = l, 

(«n<7i)i = < l-e~ s , * = 2, 

[ e 1 -^ 1 '^!, 3 < z < n, 

(H n qj) i = e- 1+( - i - 1 >D j , 1 < i < 2j — 3, j = 2, 3, . . . , m — 1, 

(66) 



(Hnqjhj-2 = 2 ' J = 2,3,...,m- 1, 



-l + e- 2 ^ + 2 S 
[Hnqjhj-i = , J = 2,3, ...,m - 1, 

(H n qjhj = CHnqjhj-2, j = 2,3, . . . ,m - 1, 
(W„ft),- = e 1 -^" 1 ^, 2j + 1 < i < n, j = 2, 3, . . . , m — 1, 
e -i+(i-i) s£)w) 1 < i < n- — 1, 
1 — e _s , i = n, 



(rl n q 
where 



/%3 



1 + e 2s - 2s 
2^ : 



Q := / e y i p j (y)dy=- 7T 1 ~> J = 2,3,... ,m - 1, 

(67) 

e(2s-l + e- 2;i ) 



e~ V( Pi(y)dy = 



2s 



D r .= e- y v J (y)dy= K — —L, j = 2, 3, . . . , m - 1, (68) 

- 1 -l + e 2s -2s 



An := / e y tp m (y)dy 



Theorem 3.1. Consider qj defined in (|19|) y>j defined in (|60l) . Let i/ie 
collocation points Xj, j = 1,2, ...,n, 6e defined in (|28[) toztft. an even number 
n > 6. TTien i/ie vectors 7-L n qj, j = — 1, 0, 1, 2, . . . , to, m = | + 1, s = ^, are 
linearly independent, where T-L n is defined in (|20[) . 

Proof. Let 

Vb := {H n q-i, H n q }, 

V s := Vj-i U {H n %}, j = l,2,...,m. 
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We prove that the elements of the sets Vj, j — 0, 1, . . . , m, are linearly indepen- 
dent. 

The elements of the set Vq are linearly independent. Indeed, W n qj 7^ Vj, 
and assuming that there exists a constant a such that 

H n q-i = aU n qo, (70) 

one gets a contradiction: consider the first and the n-th equations of (|70p . i.e., 

(H„q-i)i = Oi(H n qo)i (71) 

and 

(H n q-i) n = a(H n q )„, (72) 
respectively. It follows from (|64|) . (1651) and (|7Tj) that 

a = e 2 . (73) 

From (172 11 . ([M il . ([65 ]) and (|73 l it follows that 



- , (74) 



e 1 = e 



This is a contradiction, which proves that "H„<7__i and T-L n qo are linearly inde- 
pendent. 

Let us prove that the element of the set Vj are linearly independent, 7 = 
1,2, 3,...,m— 2. Assume that there exist constants a^, k — 1, 2, . . . , j ' + 1, such 
that 

3-1 

^™<7j = X! a k+2'H n qk- (75) 
fe=-i 

Using relations ((64|) - (l66| one can write the (2j — l)-th equation of linear system 
(|75)) as follows: 

3-1 

(H n qj)2j-i = ^2 a k+2 {H n q k ) 2] -i 

(76) 

= aie-^-V* + a 2 e- 2 +^-^ s + e 1 -^' 2 ^ £ a k+2 C k . 

k=l 

Similarly, by relations (|64[) -([66 p the (n — l)-th and n-th equations of linear 
system (I75[) can be written in the following expressions: 

4-1 

(«„<&)n-i = e- 1+2s Cj = a ie - 2+2s + a 2 e- 2s + e - 1+2s £ a k+2 C k (77) 

k=l 

and 

4-1 

CH nqj ) n = e- l+s C 3 = ai e- 2+s + a 2 e- s + er 1+s £ a k+2 C k , (78) 

fc=i 
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respectively. Multiply (1751) by e s and compare with (|77|) to conclude that a 2 = 0. 
From (|78|) with «2 = one obtains 

3-1 

at = eCj - ot k+2 C k . (79) 
fc=i 



Substitute ot\ from (|79|) and a 2 = into ([76]) and get 

(Hnqjfo-^J-W-Q'Cj. (80) 

From (pTj) and one obtains for < s < 1, j = 1, 2, 3, . . . , m — 2, the 
following relation 

e Gj- - (H n qjhj-i ^ 

(81) 

e 2s _ e -2s _ 4s S i n h(2s) - 2s V 7 

= 25 = s >0 ' 

which contradicts relation (|80"1) . This contradiction proves that the elements of 
the set Vj are linearly independent, j = 1, 2, 3, . . . , m — 2, for < s < 1. 

Let us prove that the elements of the set V m -\, are linearly independent. 
Assume that there exist constants a*;, k — 1,2, ... ,m, such that 



HnQm-l = 2J a k + 2 H n q k - (82) 
fc=-l 

Using ()64[) -(f66 |) . the (n — 3)-th equation of (|82|) can be written as follows: 

m— 3 

CH n q m -l)n-3 =e 1_4S Ari-l = a k+2 (U n q k ) n -3 

fe=-l 

m— 3 

=aie" 2+4s + a 2 e~ 4s + e~ 1+4s ^ a fc+2 C* fc + a m (H n q m - 2 ) n -3- 

k=l 

(83) 



Similarly we obtain the (n — 2)-th, (n — l)-th and n-th equations, corresponding 
to vector equation ([82]) : 



m— 3 



{U n q m -i) n -2 = one 2+3s + a 2 e 3s + e 1+3s ^ a k+2 C k + a m {U n q m -2)n-2, 



fc=i 



(84) 



m-2 



{U n q m -i)n-l = a ie - 2+2s + a 2 e- 2s + e~ 1+2s ]T a fc+2 C fc (85) 



fc=i 
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and 



CH n q m -i)n = aie 2+8 + a 2 e s + e 1+8 ^ a k+2C k , 



k=l 



respectively. Multiply ([5o]) by e s and compare with (|55|) to get 

(n n q m -i) n -i - e s {n n q m ^) n 1 - e 2s + Ase 2s - 2se 3s 



a 2 



e -2s _ I 



2s(l - e 2s ) 



(86) 



(87) 



where formula (I66[) was used. Multiplying (|86l) by e 3s , comparing with equation 
and using (j5?| . we obtain 



e'-^D^ - e ia (H n g m -i) n - ^{e^ - e 2s ) 

(H n qm-2)n-3 - e~ 1+4s C m -2 

2 + As- 2e s s + Ae 2s s - 2e 3s s + e 4s (-2 + 4s) 
-1 + e 4s - 4e 2s s ' 



(88) 



Another expression for a m is obtained by multiplying (|86[) by e 2s and comparing 
with (1541: 



(H n q m -i)n-2 - e 2s {U n qm-i)n - Me ~ e ") 

("H„g m _ 2 ) n -2 - e~ 1+3s C m -2 
2 + As- Ae s s + e 2s (-2 + 4s) 
-l + e 2s -2e s s 



(89) 



where a 2 is given in ([87]). 

In deriving formulas (|88|) and (|89| we have used the relation m = ~ + 1 and 
equation (f6"6l . Let us prove that equations ([55]) and (|59l lead to a contradiction. 
Define 



2 + As- 2e s s + Ae 2s s - 2e ds s + e 4s (-2 + 4s), 

-1 + e 4s - 4e 2s s, 

2 + As- Ae s s + e 2s (-2 + As), 

-1 + e 2s - 2e s s. 



(90) 



Then from (|88j} and (J89J) we get 

^3^2 - ri?' 4 = 0. (91) 

We have 

r 3 r2 - rir 4 = 2e s (-l + e s ) 2 s(3 + 4s + (4s - 3)e 2s - 2se s ) > for s e (0, 1). 

(92) 

The sign of the right side of equality (1921) is the same as the sign of 3 + 4s + (4s — 
3)e 2s -2se s := /3(s). Let us check that /3(s) > for s G (0, 1). One has /3(0) = 0, 
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/3'(0) = 0, P'{a) = 4-2e 2s + 8se 2s -2e s -2se s , /3" = 4e 2s + 16se 2s -4e s -2se ;i > 
0. If P"(s) > for s e (0,1) and /3(0) = 0, /3'(0) = 0, then /3(s) > for 
s e (0,1). Inequality (|52"|) contradicts relation (1§T1) which proves that H n qj, 
j = — 1,0,1,2,..., Til — 1, are linearly independent. 

Similarly, to prove that H n qj, j = — 1, 0, 1, 2, . . . , m, are linearly indepen- 
dent, we assume that there exist constants a k , k = 1, 2, . . . , m + 1, such that 



m— 1 



^ a k+2 H n q k - (93) 
Using formulas (164l) - (|66[) . one can write the {n — 5)-th equation: 

m — 1 



k=-l 

m — 4 

a ie - 2+6s + a 2 e" 6s + e - 1+6s a ]+2 C 3 

i=i 

-,1— 6s n i „1 — 6s ; 



(94) 



77- 



n ym — 3 

Similarly one obtains the (n — 4)-th, (n — 3)-th, (n — 2)-th, (n — l)-th and 
■th equations corresponding to the vector equation 

m— 4 

(H n q m )n-4 = e l -*> s D m = a ie ~ 2+5s + + e - 1+5s ^ a J+2 C 3 

l( - Hn<Zm-3)n-4 + 2 )n-4 + "m+ie S £> m _i, 

(95) 



m— 3 



/•u \ l-4s r> -2+4s , -4s . -l+4s r< 

(n n q m )n-3 = e D m = a x e + a 2 e +e > y o^C? 

i=i 

n ym— 2)71-3 +a TO+ ie s D m _i, 



(96) 



m— 3 

{U n q m ) n -2 = e l - Zs D m = a ie - 2+3s + a 2 e- 3s + e - 1+3s £ a J+2 C, 

3=1 

+ ct m ('H n q m - 2 )n- 2 + a m + i('H n q m -i) n - 2 , 

m-2 

(H n q m ) n -i = ^ 2s D m = ai e~ 2+2s + a 2 e~ 2s + e~ 1+2s £ a i+2 C,- 

3=1 

77.^771—1 jn- 1 

and 

m-2 

(%n?m)n = aie~ 2+s + a 2 e~~ s + e~ 1+s ^ otj+iCj + a m+ \(H n q m -i) n , (99) 

3=1 



(97) 



(98) 
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respectively. Here we have used the assumption n > 6. From one gets 

m-2 

«i = ('H„g m ) n e 2_s - a 2 e 2 ~ 2s - e 2J ^+2^ - a m+ \ ("H„g m -i)„e 2_s . (100) 

fc=i 

If one substitutes (|100[) into equations (|98| , (|97|) and (|96| , then one obtains the 
following relations: 

a 2 =p\ -p2Q-m+i, (101) 

0-2 = p 3 - Pi^rn - P5«m+1 (102) 

and 

«2 =P6 -P70t m -PsO^m+1, (103) 

respectively, where 

nQm — 1 )n& 

Pi '•= — ; , P2 := — 



e -2s _ 1 

e 1 3s D m — (H n q m )ne 2t (^n9m-2)«-2 — e 1+3s C m _2 

^> 3s g S g 3s 

_ (H„q m -i)n-2 - (H n q m -i)ne 2s _ e^^Dm - (H n g m ) n e 3s 

3s _ ^ 4s 

(^n9m-2)n-3 — e 1+4s C TO _2 e 1 4s _D m _i — (H n q m _i) n e 3 

P7 ■= " 



Another formula for a\ one gets from equation (|96j) : 

m— 3 

ai = e 3 ~ 8s D m - a 2 e 2 ~ Ss - e ^ ctj+iCj - a m (H n q m ^ 2 )n-3e 2 ~ ' 

3=1 



(104) 



(105) 



-a m+ ie 3 8s Z? m _i. 
Substituting (|105[) into equations (|95|) and (|94|) . yields 

"2 = P9 ~ PlO^m-l - Pll^m ~ P^m+l (106) 

and 

&2 = P9 - PU^m-l - Pl4"m ~ Pl2"m+1 , (107) 

respectively, where 

{H n q m ~3)n~A _ G 1+ ° S Cm-3 



p 9 := eD m , p 10 



e 



(^n9m-2)n-4 — 6 s {Ji n q m -2 )n-3 

Pn := 5s _ 3s , P12 := eL> m _i, 

{H n q m -3)n-5 - e- 1+6s C m -3 e 1 - 6s J D m _ 2 - e 2s (H ngm - 2 )„-3 
P13 := — jr. — . Pi* ■= e - 6s _ e - 2s ■ 

(108) 
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Let us prove that the equations (|101[) and (|106p lead to a contradiction. 
From equations (|101[) and (| 102[) we obtain 



P3 - Pi . P2 - P5 

a m = 1 

Pi Pi 

This together with (fTU5|) and ([TUT]) yield 



-own- 



P3— PI _ P6— pi 
P4 P7 

P2-P8 _ P2 — P5 
P7 P4 



Equations (fT06| . (fT07|) and (fr09)) yield 

P14-P11 (Pz-P\ , P2-P5 



_ Pu -P11 



P4 



PlO - Pl3 PlO ^ Pl3 

This together with (| 106[) imply 

' P3 - Pi , P2 - P5 



Pi 



-ot m +i 



(109) 



(110) 



(111) 



P14-P11 

"2 =P9 - PlO hPll 

PlO ~ Pl3 

where a m +i is given in (jl 10|) . Let 



T , P14-P11 . 
M :=P9- P10 hpn 



P4 



P4 



a m +l ) — Pl2Q!m+l) ( 112 ) 



PlO -P13 



P3 - Pi , P2 - P5 



P4 



/'I 



and 



c*m+i ) -Pi2a m +i (113) 



(114) 



L 2 ■= Pi - P2"m + 1- 

Then, from ([TUT]) and (fTUo) one gets 

Li-L 2 =0. (115) 

Applying formulas ([Mt - ljptj)) in (I113[) and (II 14[) and using the relation e s = 
J^jLo ff ' we obtain 

2e 4s [l + e 2s (-l + s) + s](-l + e 2s - 2e s s) 



L-i — Lo — 



(-1 + e 2s )s[3 + 4s - 2e s s + e 2s (-3 + 4s) 



2e 4 



Eoo 
7 = 



2 j I i — i 1 ■ 



J=3 0-1)! 



Eoo 
7= 



2(^-1) 



J'=2 Jl 



J + l 



(116) 



(-1 + e 2s )s[3 + 4s - 2e s s + e 2s (-3 + 4s)] 



>0, 



because e 2s > 1 for all < s < 1, 2 J '(| - j) > for all j > 3, 2(^j - 1) > for 
all j > 2 and 3+4s-2e s s+e 2s (-3+4s) > which was proved in (|52"]). Inequality 
(|116p contradicts relation (|1 15[) which proves that W n qj, j = —1,0, 1, 2, . . . , m, 
are linearly independent. 

Theorem 13. H is proved. □ 
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4 Numerical experiments 

Note that for all w,u,v £ R™ we have 

n 

WiUiVi = v t Wu, 

2 = 1 



(117) 



where t stands for transpose and 



/ tfli . . . 
w 2 



W := 







V 



w„_i 




\ 





W n J 



(118) 



Then 



DP := \\n n (f-Rh, ' 2 



"m)\\w( n )A 



n 

J2 w i n) i(f( x i) - Rh m { Xl )f + (/'(as,) - {Rh m )'{ Xl )f 



(119) 



= [Hn(f - Bh^YWHM - Rh m ) 

+ [U n (f - (Rh^YWHjf - (Rh m )'), 



where W is defined in (|118[) with ic?j = ujj , j = 1, 2, . . . ,n, defined in ([271 
The vectors R n Rh m and H n (Rh m )' are computed as follows. 
Using ([M |) - (p6]) . the vector R n Rh m can be represented by 



T-/ Rh — <? r( m ) 



(120) 



where c' m ' = 



entries: 



/ cL7 \ 

J™) 



and S m is an n x (m + 2) matrix with the following 



(5m)i,l — \Rn1-l)ii {S m )i,2 — (RnQo)i, i — 1,2, ... ,11, 

(S m )i,j = (H n qj-2)i, i = l,2,...,n, j = 3,4, . . .,m + 2. 
The vector R n {Rh m )' is computed as follows. Let 



(121) 



' '' / e v (p 3 (y)dy-e Xl / e v (pj(y)dy, i = 1, 2, . . . , n-1, j = 1, 2, . . . , m. 

(122) 
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This together with (J5DJ) yield 



i e(2i-l+ e - 2 ') . _ 1 

2/ ' 5 

2 ; e ' — — °i 
c ai e (-1+e ) 



2+e _3i -3e~ 



2/ 



= < o, 



11 



2+e 



-x< e-^-^i-l+e 2 ') 2 ■ . „".' . -. 

-e x * > »>2j + l, 

e ^ ( _i +e 2._ 20 

2 ; e > 4 ^ " •■■) 

e-'C-l+e'-Z) 



-3e~ 



i < 2j - 3 

i = 2.7 - 2 

* = 2i - 1 

* = 2.?; 



1 < i < n, 1 < j < m, 



Ji,m — 



where I = —. 

n 

Then, using (|123[) . the vector H n (Rh m )' can be rewritten as follows: 

Hn(Rh m ) = T m c 



- j 1 r ( m ) 



(123) 



(124) 



where c' m ' = 



(m) 

c o 

„( m ) 



and T m is an n x (to + 2) matrix with the following 



entries: 



{T m )i,x = -e-^ +Xi \ (T m ) i>2 = e-C 1 "**), i = l,2,...,n, 
(^m)i,j = 2- i = 1,2, ...,n, j = 3,4, . ..,m + 2. 

We consider the following examples discussed in [7]: 



(125) 



(1) f(x) = — 2 + 2cos(7r(a; + 1)) with the exact solution h(x) = — 1 + (1 + 
it 2 ) cos(7r(a; + 1)). 

(2) f(x) = -2e x_1 + - sin(7r(a; + 1)) + 2 cos(7r(izi + 1)) with the exact solution 
h(x) = i sin(vr(a; + 1)) + (1 + tt 2 ) cos(vr(a; + 1)). 

(3) f(x) = cos(^i)) + 4cos(27r(a; + 1)) - 1.5 cos( 77r( * +1) ) with the exact 
solution h(x) = |(1 + ^) cos(^il) + (2 + 8tt 2 ) cos(2vr(a; + 1)) - 0.75(1 + 
12.25tt 2 ) cos(^±^) + 1.755(a! + 1) + 2.25<5(.t - 1). 



(4) f(x) = e x + 2sin(27r(a; + 1)) with the exact solution h{x) = (1 
4tt 2 ) sin(27r(a; + 1)) + (e - 2n)S(x + 1) + 2nd(x - 1). 
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In all the above examples we have / € C 2 [— 1, 1]. Therefore, one may use the 
basis functions given in (|60[) . In each example we compute the relative pointwisc 
errors: 

\g m {U) - g(U)\ 



RPE(U).- 

maxi<i<M \g(U)\ 

where g(x) and g m (x) are defined in (j9]) and (jlip . respectively, and 

2 



*i := -1 + (i - 1)- 
The algorithm can be written as follows. 



i = l,2,...,M. 



(126) 



(127) 



Step 0. Set k = 3, n = 2k, m = f + 1, e e (0, 1) and DP > 10, where DP is 
defined in ([TT9f . 

Step 1. Construct the weights j = 1, 2, . . . , n, defined in (|27)) . 

Step 2. Construct the matrix A m+ 2 and the vector F m +2 which are defined in 

(m) 



Step 3. Solve for c 



Step 4. Compute 



( \ 

r ( m ) 



the linear algebraic system A m+ 2C = Pm+2- 



dp = \\n n (f - Rh V J 



ml\\ w (.n) 1- 



Step 5. If DP > e then set k = k + 1, n = 2k and m = § 4 
Otherwise, stop the iteration and use h m (x) = J2'j 



(128) 

1, and go to Step 1. 



fj{x) as the 

approximate solution, where (p-i(x) := 5(1 + x), ifio(x) := S(x — 1) and 
<Pj(x), j = 1,2, ... ,m, are defined in (|60p and cj m % j = —1,0,1, ... ,m, 
are obtained in Step 3. 

In all the experiments the following parameters are used: M ~ 200 and 
e = 10~ 4 , 10~ 6 , 10~ 8 . We also compute the relative error 



RE := max RPE(U 

Ki<M 



(129) 



where RPE is defined in (|126[) . Let us discuss the results of our experiments. 

Example 1. In this example the coefficients a_i and ao, given in ([S]) and ©, 

respectively, are zeros. Our experiments show, see Table 1, that the approximate 
coefficients ci^ and converge to a_i and ao, respectively, as e — > 0. Here 
to get DP < 10~ 6 , we need n = 32 collocation points distributed uniformly in 
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Table 1: Example 1 



n 


m 


e 




a_i 




m 
c -l 


a 


m 
c 




24 


13 


10" 


-4 







8.234 x 10" 5 





-5.826 x 10 


-3 


32 


17 


10" 


-6 







3.172 x 10" 5 





-1.202 x 10 


-3 


56 


29 


10" 


-8 







3.974 x 10" 6 





-6.020 x 10 


-5 


n 


m 


e 




DP 




RE 








24 


13 


10" 


-4 


8.610 x 10 


-6 


3.239 x 10" 2 








32 


17 


10" 


-6 


7.137 x 10 


-7 


1.554 x 10^2 








56 


29 


10" 


-8 


6.404 x 10 


-9 


4.337 x 10" 3 









the interval [—1,1]. Moreover, the matrix A m +2 is of the size 19 by 19 which 
is small. For e = 10~ 6 the relative error RE is of order 10~ 2 . The RPE at 
the points tj are distributed in the interval [0, 0.018) as shown in Figure 2. In 
computing the approximate solution h m at the points U, i = 1,2,..., M, one 
needs at most two out of m = 17 basis functions ifj(x). The reconstruction 
of the continuous part of the exact solution can be seen in Figure 1. One can 
see from this Figure that for e = 10~ 6 the continuous part g(x) of the exact 
solution h(x) can be recovered very well by the approximate function g m (x) at 
the points tj, j = 1,2,..., M. 




Figure 2: A reconstruction of the continuous part g{x) (above) of Example 1 
with e = 10~ 6 , and the corresponding Relative Pointwise Errors (RPE) (below) 

Example 2. This example is a modification of Example 1, where the constant 
—2 is replaced with the function — 2e x_1 + - sin(7r(:r + 1)). In this case the 
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Table 2: Example 2 



n 


m 


e 




a_i 




m 
c -l 


a 


m 
c 




24 


13 


10" 


-4 







1.632 x 10~ 4 





-1.065 x 10 


-2 


32 


17 


10" 


-6 







5.552 x 10" 5 





-2.614 x 10 


-3 


56 


29 


10" 


-8 







6.268 x 10" 6 





-1.954 x 10 


-4 


n 


m 


e 




DP 




RE 








24 


13 


10" 


-4 


8.964 x 10 


-6 


3.871 x 10" 2 








32 


17 


10" 


-6 


7.588 x 10 


-7 


1.821 x 10^2 








56 


29 


10" 


-8 


6.947 x 10 


-9 


4.869 x 10" 3 









coefficients a_i and ao are also zeros. The results can be seen in Table 2. As in 
Example 1, both approximate coefficients and c "'' ) converge to as e — >• 0. 
The number of collocation points at each case is equal to the number of col- 
location points obtained in Example 1. Also the RPE at each observed point 
is in the interval [0, 0.02). One can see from Figure 3 that the continuous part 
g(x) of the exact solution h(x) can be well approximated by the approximate 
function g rn (x) with e = 10~ 6 and RE = 0(1Q- 2 ). 



jpc-azriiilkn i;i m > ve e/.if.l ijj 




Figure 3: A reconstruction of the continuous part g(x) (above) of Example 2 
with e = 10 -6 , and the corresponding Relative Pointwise Errors (RPE) (below) 

Example 3. In this example the coefficients of the distributional parts a_i and 
ao are not zeros. The function / is oscillating more than the functions / given 
in Examples 1 and 2, and the number of collocation points is larger than in the 
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Table 3: Example 3 



n 


m 


e 


a_i 


m 
C -l 


a 




80 


41 


10" 4 


1.750 


1.750 


2.250 


2.236 


128 


65 


10~ 6 


1.750 


1.750 


2.250 


2.249 


232 


117 


10" 8 


1.750 


1.750 


2.250 


2.250 


n 


m 


e 


DP 


RE 






80 


41 


10" 4 


4.635 x 10" 5 


2.282 x 10" 2 






128 


65 


10" 6 


9.739 x 10~ 7 


7.671 x 10~ 3 






232 


117 


10" 8 


7.804 x 10" 9 


2.163 x 10" 3 







previous two examples, as shown in Table 3. In this table one can see that the 
approximate coefficients and converge to a_i and a 0l respectively. The 
continuous part of the exact solution can be approximated by the approximate 
function g m (x) very well with e = 10~ 6 and RE = O(10 -3 ) as shown in Figure 
4. In the same Figure one can see that the RPE at each observed point is 
distributed in the interval [0,8 x 10~ 3 ). 



npp-azi-ialb-i !a M| :i vs ig) 




Figure 4: A reconstruction of the continuous part g(x) (above) of Example 3 
with e = 10~ 6 , and the corresponding Relative Pointwise Errors (RPE) (below) 



Example 4. Here we give another example of the exact solution h having 
non-zero coefficients a_i and oo- In this example the function / is oscillating 
less than the / in Example 3, but more than the / in examples 1 and 2. As 
shown in Table 4 the number of collocation points n is smaller than the the 
number of collocation points given in Example 3. In this example the exact 
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coefficients a_i and ao are obtained at the error level e = 10~ 8 which is shown 
in Table 4. Figure 5 shows that at the level e = 10~ 6 we have obtained a good 
approximation of the continuous part g(x) of the exact solution h(x). Here the 
relative error RE is of order O(10~ 2 ). 



Table 4: Example 4 



n 


m 


e 




a_i 




Am) 
c -l 


ao 


c 


40 


21 


10" 


4 


-3.565 




-3.564 


6.283 


6.234 


72 


37 


10" 


6 


-3.565 




-3.565 


6.283 


6.279 


128 


65 


io- 


8 


-3.565 




-3.565 


6.283 


6.283 


n 


777 


e 




DP 




RE 






40 


21 


10" 


4 


8.775 x 10 


-5 


3.574 x 10" 2 






72 


37 


10" 


6 


6.651 x 10 


-7 


1.029 x 10- 2 






128 


65 


10" 


8 


6.147 x 10 


-9 


3.199 x 10" 3 








Figure 5: A reconstruction of the continuous part g(x) (above) of Example 4 
with e = 10~ 6 , and the corresponding Relative Pointwise Errors (RPE) (below) 
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